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HIGH-ORDER NONLINEAR 

SCHRODINGER EQUATION FOR THE ENVELOPE 
OF SLOWLY MODULATED GRAVITY WAVES 
ON THE SURFACE OF FINITE-DEPTH FLUID 
AND ITS QUASI-SOLITON SOLUTIONS 


We consider the high-order nonlinear Schrodinger equation derived earlier by Sedletsky 
[Ukr. J. Phys. 48(1), 82 (2003)] for the first-harmonic envelope of slowly modulated gravity 
waves on the surface of finite-depth irrotational, inviscid, and incompressible fluid with flat 
bottom. This equation takes into account the third-order dispersion and cubic nonlinear disper¬ 
sive terms. We rewrite this equation in dimensionless form featuring only one dimensionless 
parameter kh, where k is the carrier wavenumber and h is the undisturbed fluid depth. We show 
that one-soliton solutions of the classical nonlinear Schrodinger equation are transformed into 
quasi-sohton solutions with slowly varying amplitude when the high-order terms are taken into 
consideration. These quasi-soliton solutions represent the secondary modulations of gravity 
waves. 
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1. Introduction 


( 1 . 1 ) 


The nonlinear Schrodinger equation (NLSE) 

A r = — Qj\A x — ia 2 A xx -\- iao, o, oA|A| 2 

arises in describing nonlinear waves in various phys¬ 
ical contexts, such as nonlinear optics j(54| . plasma 
physics [dill] , nanosized electronics jl2l , ferromagnet¬ 
ics 0 , Bose-Einstein condensates [73) , and hydrody¬ 
namics El S3, m,Si I- Here, x I s the direction of wave 
propagation, r is time, A(\, t) is the complex first- 
harmonic envelope of the carrier wave, and the sub¬ 
scripts next to A denote the partial derivatives. NLSE 
takes into account the second-order dispersion (term 
with A xx ) and the phase self-modulation (term with 
H|H| 2 ). The coefficients a±, a 2 , and ao, 0,0 take vari¬ 
ous values depending on the particular physical con¬ 
text under consideration. 

In the general context of weakly nonlinear disper¬ 
sive waves, this equation was first discussed by Ben- 
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ney and Newell [5] - In the case of gravity waves prop¬ 
agating on the surface of infinite-depth irrotational, 
inviscid, and incompressible fluid, NLSE was first de¬ 
rived by Zakharov |68j using the Hamiltonian formal¬ 
ism and then by Yuen and Lake [66j using the aver¬ 
aged Lagrangian method. The finite-depth NLSE of 
form (11.11) was first derived by Hasimoto and Ono |3Q| 
using the multiple-scale method and then by Stiassnie 
and Shemer [52] from Zakharov’s integral equations. 
Noteworthy is also the recent paper by Thomas et 
al. IHj who derived the finite-depth NLSE for water 
waves on finite depth with constant vorticity. 

Under certain relationship between the parameters, 
when 

02(10,0,0 < 0 , ( 1 - 2 ) 

NLSE admits exact solutions in the form of solitons, 
which exist due to the balance of dispersion and non¬ 
linearity and propagate without changing their shape 
and keeping their energy 0. In this case, the uni¬ 
form carrier wave is unstable with respect to long¬ 
wave modulations allowing for the formation of en- 
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velope solitons. This type of instability is known as 
the modulational or Benjamin-Feir instability Q (it 
was discovered for the first time in optics by Bespalov 
and Talanov 0 ) . In the case of surface gravity waves, 
condition holds at kh > 1.363, k being the car¬ 
rier wavenumber and h being the undisturbed fluid 
depth. In addition to theoretical predictions, enve¬ 
lope solitons were observed in numerous experiments 
performed in water tanks 0, 44, 5(| 5j], 0,Ek],EiB,E0 ■ 
At the bifurcation point a o,o,o = 0 (kh ~ 1.363), 
when the modulational instability changes to stabil¬ 
ity, NLSE of form (11.11) is not sufficient to describe 
the wavetrain evolution since the leading nonlinear 
term vanishes. In this case, high-order nonlinear and 
nonlinear-dispersive terms should be taken into ac¬ 
count. In the case of infinite depth, such a high-order 
NLSE (HONLSE) was first derived by Dysthe [H]. It 
includes the third-order dispersion (A xxx ) and cubic 
nonlinear dispersive terms (|A| 2 A X , A 2 A *, asterisk 
denotes the complex conjugate) as well as an addi¬ 
tional nonlinear dispersive term describing the input 
of the wave-induced mean flow (some of these terms 
were introduced earlier by Roskes f4(| without tak¬ 
ing into consideration the induced mean flow). This 
equation is usually referred to as the fourth-order 
HONLSE to emphasize the contrast with the third- 
order NLSE. Janssen [34] re-derived Dysthe’s equa¬ 
tion and corrected the sign at one of the nonlinear dis¬ 
persive terms. Hogan [31! ] followed the earlier work by 
Stiassnie [56] to derive the similar equation for deep¬ 
water gravity-capillary waves with _surface tension 
taken into account. Selezov et al. 0] extended the 
HONLSE derived by Hogan to the case of nonlinear 
wavetrain propagation on the interface of two semi¬ 
infinite fluids without taking into account the induced 
mean flow. Worthy of mention is also the paper by 
Lukomsky [HJ who derived Dysthe’s equation in a dif¬ 
ferent way. Later, Trulsen and Dysthe §32] extended 
the equation derived by Dysthe to broader bandwidth 
by including the fourth- and fifth-order linear dis¬ 
persion. Debsarma and Das 0 derived a yet more 
general HONLSE that is one order higher than the 
equation derived by Trulsen and Dysthe. Gramstad 
and Trulsen [0 derived a set of two coupled fourth- 
order HONLSEs capable of describing two interact¬ 
ing wave systems separated in wavelengths or di¬ 
rections of propagation. Zakharov and Dyachenko 
made a conformal mapping of the fluid 
domain to the lower half-plane to derive a counter¬ 


part of Dysthe’s equation in new canonical variables 
(the so-called compact Dyachenko-Zakharov equa¬ 
tion f2lj|2S|). 

Original Dysthe’s equation was written for the first- 
harmonic envelope of velocity potential rather than of 
surface profile. In the case of standard NLSE, this dif¬ 
ference in not essential because in that order the first- 
harmonic amplitudes of the velocity potential and 
surface displacement differ by a dimensional factor 
only, which is not true anymore in the HONLSE case, 
as discussed by Hogan [33] ■ Keeping this in mind, 
Trulsen et al. [63] rewrote Dysthe’s equation in terms 
of the first-harmonic envelope of surface profile while 
taking into account the linear dispersion to an arbi¬ 
trary order. 

In the case of finite depth, the effect of induced 
mean flow manifests itself in the third order, so that 
the NLSE is generally coupled to the equation for 
the induced mean flow 0. However, Davey and Stew- 
artson [13] showed that these coupled equations are 
equivalent to the single NLSE derived by Hasimoto 
and Ono (30j |. On the other hand, such an equiva¬ 
lence is not preserved for high-order equations. The 
first attempt to derive a HONLSE in the case of fi¬ 
nite depth was made by Johnson j35j|, but only for 
kh s=s 1.363, when the cubic NLSE term vanishes. The 
similar attempt was made by Kakutani and Michihiro 
0 (see also a more formal derivation made later by 
Parkes ®). A general fourth-order HONLSE for the 
first-harmonic envelope of surfaceprofile was derived 
by Brinch-Nielsen and Jonsson [8j in coupling with 
the integral equation for the_wave-induced mean flow. 
Gramstad and Trulsen 00 derived a fourth-order 
HONLSE in terms of canonical variables that pre¬ 
serves the Hamiltonian structure of the surface wave 
problem. 

Sedletsky 0,0] used the multiple-scale technique 
to derive a single fourth-order HONLSE for the first- 
harmonic envelope of surface profile by introducing 
an additional power expansion of the induced mean 
flow. This equation is the direct counterpart of Dys¬ 
the’s equation written in terms of the first-harmonic 
envelope of surface profile [63] but for the case of fi¬ 
nite depth. Slunyaev [0] confirmed the results ob¬ 
tained in [0] and extended them to the fifth order. 
Grimshaw and Annenkov 0] considered a HONLSE 
for water wave packets over variable depth. 

The deep-water HONLSE in the form of Dysthe’s 
equation was extensively used in numerical simula- 
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tions of wave evolution Mi 
However, no such modeling has been performed in 
the case of finite depth because of the complexity of 
equations as compared to the deep-water limit. The 
equation derived in [HHl. [Till can be used as a good 
starting point for the simulations of wave envelope 
evolution on finite depth. The aim of this paper is (i) 
to rewrite this equation in dimensionless form suit¬ 
able for numerical integration and (ii) to observe the 
evolution of NLSE solitons taken as initial waveforms 
in the case when the HONLSE terms are taken into 
consideration for several values of intermediate depth. 

This paper is organized as follows. In Section [21 we 
write down the fully nonlinear equations of hydro¬ 
dynamics used as the starting point in this study. In 
Section [3] we formulate the constraints at which the 
fully nonlinear equations can be reduced to HONLSE. 
Then we briefly outline the multiple-scale technique 
used to derive this equation, which is presented in 
Section 0] Next, we introduce dimensionless coordi¬ 
nate, time, and amplitude to go over to the dimen¬ 
sionless HONLSE. As a result, only one dimension¬ 
less parameter kh appears in the equation. The final 
step is to pass to the reference frame moving with 
the group speed of the carrier wave. In Section [5l we 
present the results of numerical simulations and com¬ 
pare the NLSE and HONLSE solutions. Conclusions 
are made in Section [6] 


II m, El [Hi [Hi |. 


2. Problem Formulation 


We consider the dynamics of potential two-dimen¬ 
sional waves on the surface of irrotational, inviscid, 
and incompressible fluid under the influence of grav¬ 
ity. Waves are assumed to propagate along the hori¬ 
zontal cc-axis, and the direction of the vertical y-axis 
is selected opposite to the gravity force. The fluid is 
assumed to be bounded by a solid flat bed y = —h 
at the bottom and a free surface y = r)(x, t) at the 
top (Fig. [T]). The atmospheric pressure is assumed to 
be constant on the free surface. Then the evolution 
of waves and associated fluid flows is governed by the 
following system of equations 24,58]: 


+ & yy = 0, —oo < x < oo, —h < y < r]{x , t); 

( 2 . 1 ) 

+ + +gv = 0, y = r)(x,t); (2.2) 

r]t- ® y + rix®x = 0, y = rj(x, f); (2.3) 



Fig. 1. Sketch of the physical domain occupied by an ideal 
incompressible fluid of finite depth 


$ y = 0, y=-h-, (2.4) 

where <J>(a;, y , t ) is the velocity potential (the velocity 
is equal to Vd*), g is the acceleration due to gravity, 
t is time. Here, (12.11) is the Laplace equation in the 
fluid domain, (12.21) is the dynamical boundary con¬ 
dition (the so-called Bernoulli or Cauchy-Lagrange 
integral), (12.31) and (12.41) are the kinematic boundary 
conditions (no fluid crosses the free surface and the 
bottom), the indices x, y, and t designate the partial 
derivatives over the corresponding variables. The po¬ 
sition of the zero level y = 0 is selected such that the 
Bernoulli constant (the right-hand side of Eq. (12.21) 1 
is equal to zero. 

Consider a modulated wavetrain with carrier fre¬ 
quency uj and wavenumber k. In this case, a solution 
to Eqs. (HOT) (l2~4l) can be looked for in the form of 
Fourier series with variable coefficients: 

/$(x, y, t)\ _ \ ' Vi t)\ ±n(kx—ujt) 

\ »?(*. *) / n i" oc V Vnix, t) ) ’ (2.5) 

V-n = Vn, $-n = $0, 

where * stands for complex conjugate (here, we as¬ 
sume the carrier wave to be symmetric), the functions 
c f > (x, y, t ) and r)(x, t ) are assumed to be real by defi¬ 
nition. Substituting (12.51) in (12.11) (12.41) and equating 
the coefficients at the like powers of exp(i(fcx — wf)), 
one can obtain a system of nonlinear partial dif¬ 
ferential equations for the functions $,, ( 2 ;, y , t) and 
r] n (x, t). Linearization of these equations at n = 1 
gives the dispersion relation for gravity waves: 

w 2 = gk tanh(fch). (2-6) 
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3. Slowly Modulated Quasi-Harmonic 
Wavetrains and Multiple-Scale Expansions 

Generally the system of equations for & n (x, y, t) 
and q n {x, t) is by no means more simple than orig¬ 
inal equations. It can be simplified when solutions 
are looked for in a class of functions with nar¬ 
row spectrum, |Afc| <C k. In this case, the problem 
has a formal small parameter y ~ |Afc|/A; (quasi¬ 
monochromaticity condition), with <F„(;r, y, t) and 
rj n (x, t ) being slow functions of x and t. Accordingly, 
the wave motion can be classified into slow one and 
fast one by introducing different time scales and dif¬ 
ferent spatial scales: 


T i Tlj. x/" _ Tl _ 

n = fd t) X-n — fd X. 


(3.1) 


The derivatives with respect to time and coordinate 
are expanded into the following series: 


d_ 

dt 


OO 


71=0 


d 

dT n ’ 


d_ 

dx 


OO 


71=0 


d 

dx n : 


(3.2) 


the times T n and coordinates X n being assumed to 
be independent variables. 

When there are no resonances between higher har¬ 
monics, the amplitudes of Fourier coefficients de¬ 
crease with increasing number (quasi-harmonicity 
condition): 


ij n ~ e n A, n ^ 1, r]o ~ e 2 A, e < 1, (3.3) 


where 

Vi = ^ £ A(x, t). (3.4) 


The parameter e can be regarded as a formal small 
parameter related to the smallness of wave amplitude 
as compared to the carrier wavelength A = ^. In this 
case, the unknown functions $„(:e, y, t ) and i] n (x, t ) 
can be expanded into power series in the formal pa¬ 
rameter e: 



Multiple-scale expansions (13.21) and m allow the 
functions $ n (a;, y, t) and r] n (x, t) to be expressed in 
terms of the first harmonic envelope A(x, t), as de¬ 
scribed in detail in [43] • Note that in the procedure 
described in 47| it is essential to set e = y. 


In practice, the quasi-harmonicity condition can be 
written as 


Nil < 1, 


(3.6) 


and the condition of slow modulation (quasi-mono¬ 
chromaticity) can be formalized as 


A x 

kA 


<C 1, 


(3.7) 


which follows from differentiating the function 
A(x, t ) exp(i (kx — ut)) with respect to x. With these 
conditions satisfied, the original system of equations 
(12.11) (12.41) can be reduced to one evolution equation 
for the first harmonic envelope A{x, t) with the use 
of small-amplitude expansions (13.21) and (13.51) . 


4. High-Order 

Nonlinear Schrodinger Equation 


4-1. Equation derived by Sedletsky 

Sedletsky 0,0 used the above-described multiple- 
scale procedure to derive the following HONLSE for 
the first-harmonic envelope A(x, t) (Eq. (68) in [47 


if 


dA 


V„ 


dA 


\d(et) y d(ex) 


+ 


+e + “*’ 2 ® 1 ' 412 ' 1 ) + 

1 A ■ ukQ„{Af 3A 


1£ —-W 


-{- (xkQ^2 A 


6 d(ex) 3 
2 dA 


d{sx) 


d(ex] 


= 0 [m/s]. 


(4.1) 


As compared to the standard NLSE, this equation 
takes into account additional nonlinear and disper¬ 
sive terms of order 0(g 2 ). Equation (14.11) was later 
re-derived by Slunyaev [521 ], who confirmed the sym¬ 
bolic computations presented in 0,0 and extended 
them to the 0(e 3 ) order. Here, we restrict our atten¬ 
tion to the original equation m- The parameters of 
this equation are given by 


w = (gka) 1/2 , 


a = tanh(fch), 


du_ w / 2 kh \ 

dk 9 2k \ sinh(2 kh)J 


to 

2k 



(4.2a) 


(4.2b) 
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d 2 


UJ 


UJ 


UJ = 


dk 2 4k 2 a 2 

— 2a ( a 2 — l) kh — a 2 ) 


((a 2 - 1) (3a 2 + 1) k 2 h 2 - 


(4.2c) 


UJ = 


<9 3 w 
dk 3 

.4 


U 


8 fc 3 a 3 


((a 2 - 1) x 


x (15a 4 — 2a 2 + 3) k 3 h 3 — 3a (a 2 — l)(3a 2 + l) fc 2 /i 2 — 
— 3a 2 (a 2 — 1 ) kh — 3a 3 ), (4.2d) 

93 = ~ ^ “ 10a2 + 9 ) k2h ‘ 2 + 

+ 2a (3a 6 - 23a 4 + 13a 2 - 9) kh - 


- a 2 (7a 4 - 38a 2 - 9)), (4.2e) 

Q “ = sdv (<- 2 - ‘) s x 

x (3a 6 - 20a 4 - 21a 2 + 54) k 5 h 5 - 

-a (a 2 - l ) 3 (lla 8 — 99a 6 - 61a 4 + 7a 2 + 270 )k 4 h 4 + 

+ 2a 2 (a 2 - 1) (7a 10 - 58a 8 + 38a 6 + 52a 4 - 

- 181a 2 + 270 )k 3 h 3 - 2a 3 (3a 10 + 18a 8 - 146a 6 - 
-172a 4 +183a 2 -270)fc 2 h 2 -a 4 (a 8 -109a 6 +517a 4 + 

+ 217a 2 +270)fch+a 5 (a 6 - 40a 4 + 193a 2 + 54)) +A, 
i (4.2f) 

x (3a 6 + 7a 4 - lla 2 + 9) k 5 h 5 + 

+ a (a 2 — l ) 3 (lla 8 — 48a 6 + 66 a 4 + 8 a 2 + 27) fc 4 h 4 — 

- 2a 2 (a 2 - 1) (7a 10 - 79a 8 + 282a 6 - 

- 154a 4 - a 2 + 9 )k 3 h 3 + 2a 3 (3a 10 - 63a 8 + 314a 6 - 

- 218a 4 + 19a 2 + 9 )k 2 h 2 + a 4 (a 8 + 20a 6 - 158a 4 - 

- 28a 2 —27)/c/i—a 5 (a 6 - 7a 4 + 7a 2 - 9))-A, (4.2g) 

v = (a 2 - l ) 2 k 2 h 2 - 2a (a 2 + l) kh + a 2 . (4.2h) 


The quantity V g is the wave group speed. The param¬ 
eter A is the correction introduced by Slunyaev [52f 
to the coefficients derived in [13, SH- This correction 
is negligible at kh > 1 (see Appendix [A}, and we 
ignore it by keeping A = 0. 


The free-surface displacement is expressed in terms 
of A as 


rj = £ 2 Vo + e Re(Ae 1(?cx -^ ) ) + 

+ e 2 2 Re(r ?2 e 21 ^-^) + 0(e 3 ), 


(4.3) 


where Re{-} stands for the real part of a complex¬ 
valued function. Here, rjo and 772 are defined as 


ho = 


„ + 2 (i-^)a t|J| , 


m = 


3-a 2 
8 a 3 


kA 2 . 


(4.4a) 

(4.4b) 


The corresponding velocity potential is written as 
$ = e$o + e 2 Re($i e 1 ( fe3 '-“ t )) + 

+ £ 2 2 Re ($ 2 e 2i ( fca: - wt )) + 0( e 3 ), (4.5) 

where 


$1 = 


d A( h , + ^ ] _ iA 

UJ 


~(y + h) 


2 ka \\dx 

dA sinh (k(y + h)) 


cosh(fc(y + h)) 


cosh (kh) 


dx cosh (kh) 


$ 2 = 3iw 


(a 4 - 1 ) cosh( 2 k(y + h)) a2 
16a 4 cosh( 2 fch) 


A 2 . 


(4.6a) 


(4.6b) 


The term $0 describes the wave-induced mean flow 
and is expressed implicitly in terms of its derivatives 


dx 


eAlilAf 

2av 


+ ie 


W72 
8 a 2 v 2 



d$o _ d4>o 
dt ~ 9 dx' 



(4.7a) 

(4.7b) 


where 


71 = (a 2 — 1 ) 2 kh — a(a 2 — 5), (4.8a) 

72 = (a 2 — l) 5 fc 4 h 4 + 4a(a 2 — l)'(l3a 2 + 3 )k 3 h 3 — 
- 2a 2 (a 2 - 1) (3a 4 + 32a 2 - 3) k 2 h 2 + 


+ 4a 3 (2a 4 — a 2 — 5) kh — 3 a 4 (a 2 — 5). (4.8b) 

Functions (14.311 and (14.511 define an approximate so¬ 
lution to the original system of equations EH) EH 
in terms of the first-harmonic envelope A, which is 
found from Eq. EH- 
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f.2. Dimensionless form 

Introduce the following dimensionless time, coordi¬ 
nate, and amplitude: 

r = (3t, x = kx, u = a~ 1 eA, (4.9) 


a and (3 being the parameters to be determined. The 
relationship between the old and new derivatives is 

d d d d 
fa =k d^’ dt = (j d t' 

Then Eq. (14.11) is transformed to 


ia (/ 3u t + kV g u x )+a ^-u/'fc 2 u xx + wfc 2 < 73 |a| 2 |u| 2 uj + 
+ i a (—^ UJ '"k 3 u X xx + tok 2 Q4i\a\ 2 u x \u\ 2 + 

+ ojk 2 Q 4 2 \a\ 2 u 2 u 1 x j = 0 [m/s]. 

Here, the indices \ and r designate the partial deriva¬ 
tives with respect to the corresponding variables. Ta¬ 
king into account that ui" < 0 at all h > 0 (Fig. [2]), 
divide this equation by u>"k 2 a so that 


ui"k 2 


V g 

——? 

(jj"k 


+ 1 ^ U XX + Q3\U\ U ) + 


. ( 1 U)"'k W . 2 2 

1 \ 7; ,// U xxx + -^\a\ Qiiu x \u\ + 

\ 0 CO CO 


CO 


a 


l Q42U 2 U*^j 


= 0 


and select the values of a and f3 as 


\ot\ 


— — > 0, /3 = -J'k 2 > 0. 

ui 


(4.10) 



Thus, Eq. (14.11) takes the dimensionless form 


V g \ 1 

1 - -J7^ u xj - ? u xx + q 3 \ u \-u + 

1 u}"'k 

1 1 G~uN~ Uxxx + Q ilU x\ u \ 2 + Qa 2 U 2 u* 
or, equivalently, 


= 0 


i [u T + a\u x ) — a 2 u xx + ao,o,oM 2 u + i^—d 3 U xxx + 
+ oi,o,o'OxM 2 + a o,o, i« 2 u x j =0, (4.11a) 


which finally yields 

u T = —a\u x — ici 2 U xx + iao, o, o|w| 2 u+ 

+ (a 3 u xxx - ay o, ou x \u\ 2 - a 0 ,o. \u 2 u*^j, (4.11b) 

where we used the unified notation introduced by 
Lukomsky and Gandzha [Q. Here, the coefficients 


a\ = ——? 7 T = — — (cr 2 + a (l - a 2 ) kh] > 0, 
u>"k v\ / 

1 

a 2 --, 

1 u m k 

a 3 = - — = 

6 oj 

= —-— ((a 2 — l) (15a 4 — 2 o 2 + 3) k 3 h 3 — ^ 

— 3cJ (er 2 — l) (3a 2 + l) k 2 h 2 — 

— 3a 2 (cr 2 — l) kh — 3a 3 ), 

^ 0 , 0,0 = 93) Ol, 0,0 = Qih Oo, 0 ,l = Q 42 , 

v=(cr 2 - 1) (3a 2 + 1) k 2 h 2 - 2a (a 2 - 1 ) kh - cr 2 


are all real and depend on one dimensionless pa¬ 
rameter kh. Their behavior as functions of kh is 
shown in Fig. [3] It can be seen that Eq. (14.11) is valid 
at kh > 1, where the coefficients ao,o,o> oyo.Oi and 
do, 0,1 do not diverge. At smaller depths, the Korte- 
weg-de Vries equation and its generalizations @,0] 
should be used. On the other hand, at large kh , 
the infinite-depth limit (Dysthe’s equation) should be 
used. Indeed, the following asymptotics are easily ob¬ 
tained at kh — > 00 : 

1 1 3 i 

03 — -j^> oq, 0,0 — — 2 > a i,o,o — ao,o,i — j- ( 4 - 13 ) 
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Fig. 3. Normalized coefficients of HONLSE as functions of kh 


They coincide with the corresponding coefficients of 
Dysthe’s equation [63|, except for the term includ¬ 
ing the wave-induced mean flow, which cannot be 
explicitly reconstructed from Eq. (14.111) because of 
the additional power expansion of the wave-induced 
mean flow made to derive Eq. (ED- However, this 
term can be reconstructed from the equations gene¬ 
rating Eq. ED, at the stage when the wave-induced 
mean flow has not been excluded from the equation 
for A yet [47]. Taking into account these constraints, 
we will restrict our attention to the following range 
of intermediate depths: 

1 < kh < 5. (4.14) 

4-3. Moving reference frame 

Equation (14.111) can be rewritten in the form with¬ 
out the u x term. To this end, let us proceed to the 
reference frame moving with speed ai (dimensionless 
group speed): 

£ = X ~ ai t, T = t. (4.15) 

The relationship between the derivatives in new and 
old variables is given by the formulas 

d _di d dT d _ d 
d^~ d^d£ + d^df~W 
d _ d( d dT d _ d d 
d t ~ ~(h~di + IfrdT ~ ~ ai ~dl + dT’ 


+ ( a 3 u |it| 2 - a 0 ,o, im 2 m|J- (4-16) 

This is our target equation for numerical simula¬ 
tions. It possesses the integral of motion 


OO 

t o( t ) = j K£, T)| 2 d£ = const, 

— OO 


(4.17) 


which expresses the conservation of wave action. The 
derivation of this conservation law is given in Ap- 
pendix[Cj It allows one to trace the relative numerical 
error of simulations: 


Er(/ 0 ) = 


|/o(t)-/q(0)| 

Io(0) 


(4.18) 


Of particular interest is to reveal any relationship 
of Eq. (14.1611 to other HONLSEs derived in differ¬ 
ent contexts. In Appendix m we consider one such 
equation (the Sasa-Satsuma equation) and prove that 
Eq. (14.161) cannot be reduced to it at any kh. 


4-4 ■ Dimensionless free surface 
displacement and velocity potential 

The dimensionless free surface displacement is ex¬ 
pressed in terms of u as follows: 


so that 

u T = —i a 2 ii££ + ia 0i o,oM 2 u-|- 


£ = ky = ao|rt| 2 + u\ Re(w e lS ) + 2«2 Re(n 2 e 2lS ), 

(4.19) 
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Fig. 4- Ratio between the dimensionless phase and group 
speeds as a function of kh 


cosh(z + kh) 
cosh {kh) 


(z + kh) 


sinh (2 + kh) 
cosh (kh) 



3i(er 4 — l) cosh(2 (z + kh)) 2 
16cr 4 cosh( 2 kh) U ’ 

where 2 = ky is the dimensionless vertical coordi¬ 
nate. The quasi-harmonicity condition is written as 

< 1, (4.24) 

\/c 

and the quasi-monochromaticity condition is 


cr + 2 (l — a 2 ) kh 1 3 — a 2 

ao ; c^i 0.2 1 

cv yjc oca A 

where 

9 = kx — cut = x ~ ct = t; + (ai — c) r (4.20) 

is the wave phase and 


< 1 . 


(4.25) 


Finally, the original equations of hydrodynamics can 
be written in the following dimensionless form: 

iptf + <-Pzz = 0 , -00 < £ < 00 , -kh < 2 < C(£, t); 

(4.26) 


1 _ 4 cr 2 

\a\ 2 k 2 v 


(4.21) 


is the dimensionless phase speed. Figure 2] shows the 
ratio of the phase speed c to the group speed «i as a 
function of kh. This ratio is equal to unity at kh —> 0, 
and it is twice as large at kh — > 00 , in full conformity 
with the classical water wave theory (i3 |. The wave 
envelope is written as 


Envelope ~ a l M + ( a 0 + 2 «2 ) W\ 2 ■ (4.22) 

The corresponding dimensionless velocity potential 
is expressed as 


‘fir + ^ (^C + ^z) + — C — Oj 2 — C(C; r ); (4.27) 

(t-‘Pz+Q‘P£ = o, z = C(£, t); (4.28) 

ip z = 0, 2 = — kh. (4.29) 

5. Numerical Simulations 

In this section, we adopt the split-step Fourier (SSF) 
technique described in Appendix [D] to compute solu¬ 
tions to HONLSE 114.1611 . To test the accuracy of our 
numerical scheme, we start from classical NLSE ED 
written in terms of the coordinate £. At a 0 , 0,0 < 0 
(kh > 1.363), it has an exact one-soliton solution |69j: 


<p = —= Vo + 2 Re(^i e 1 ®) + 
w" 

+ 2 Re((p 2 e 2ie ), 

^ = 2^ |u| +g^K -««£) 


(4.23) 


w(C, r) = 

n = (k 2 


u 0 exp (i«£ — iflr) 

cosh(A'(£ - £0 - Vt)) ’ 

— k 2 ) a 2 , V = —2 na 2 , = |u 0 | 


1^0 £ C, Av, ^0 ^ K- 


(5.1) 


/ 


Qp, 0 , 0 

2a 2 ’ 


(<A)) t = -ai (ipo)£ , 



(cr 2 + l)fc/i + a \ 

- 2a - U T 


Here, V is the soliton speed, Uq is the complex am¬ 
plitude, k and Q are the soliton’s wavenumber and 
frequency, and £0 is the soliton’s initial position. The 
amplitude uq and wavenumber k should be selected 
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such that the quasi-harmonicity and quasi-monochro¬ 
maticity conditions (14.241) . (14.2511 hold true. In prac¬ 
tice, these conditions mean that the soliton amplitude 
and wavenumber should be small: 

|uo| Cl, K«l. 

In this study, we restrict our attention by the follow¬ 
ing choice of parameters: 

u 0 = 0.1, k = -K (=> ft = 0), fo = 0. (5.2) 

Figures ED and ED demonstrate that constraints (14.24[) 
and (14.2511 are readily satisfied in this case. Note that 
at k < 0 we have V > 0. In this case, solitons move 
from left to right with speed exceeding the carrier 
group speed. 

Figure [7] shows a soliton computed for kh = 3 using 
analytical formula (15.11) for the initial moment r = 0 
and moment r = 10000. The same soliton was taken 
as the initial condition for the simulation with the 
SSF technique. The deviation from the exact solution 
is seen to be negligible. Indeed, the numerical error 
estimated with formula (14.181) is 

S (2) |r=ioooo : Er(J 0 ) = 1.3 x 10~ 10 %, 

^rmsf^exact, ^comp) = 1-0 X 10 %, 

S (4) |r=ioooo : Er(J 0 ) = 2.3 x 10" 10 %, 

^rmsf^exact, ^comp) ~ 3.5 X 10 %, 

where S^ and designate the order of the SSF 
technique adopted for calculation (see Appendix [Pi) 
and 


A rms (u, g){r) 


IZ (K£, t)| - \g(£, r)|) dC 


is the relative r.m.s. deviation between two func¬ 
tions. Thus, our numerical scheme reproduces the ex¬ 
act one-soliton solution to NLSE with high accuracy. 

Figure [ 8 ] shows the evolution of the same one- 
soliton waveform taken as the initial condition in 
HONLSE (14.161) . As compared to the NLSE case, the 
wave amplitude is smaller, the pulse width is larger, 
and the wave speed is higher. The wave amplitude 
does not remain constant and exhibits slow oscilla¬ 
tions that can be interpreted as the secondary mod¬ 
ulation of the carrier wave. The amplitude of these 



Fig. 5. Testing the quasi-harmonicity condition (14.24ft for 
kh = 3 



Fig. 6. Testing the quasi-monochromaticity condition Il4.25ft 
for kh = 3 


oscillations decreases with time (Fig.ED). Such a so¬ 
lution does not fall under the definition of soliton 
because it does not preserve the constant amplitude 
and shape during the evolution. On the other hand, 
it moves with nearly the constant speed (Fig. [TUI) and 
still possesses the unique property of solitons to exist 
over long periods of time without breaking. In view 
of this unique property, we call such solutions quasi- 
solitons. The term quasi-soliton was introduced ear¬ 
lier by Zakharov and Kuznetsov 0 , but in somewhat 
different context; then Karpman et al. [.38l| and Slun- 
yaev [53| used it in the same context as in the present 
study. 

Such a behavior of NLSE solitons in the HONLSE 
case was first described by Akylas J3] in the context 
of asymptotic modeling and numerical simulations of 
Dysthe’s equation in the infinite-depth limit. Growth 
in the soliton speed corresponds to the well-known 
carrier frequency downshift observed in deep-water 
experiments by Su [ 6 (j and in simulations of Dysthe’s 
equation by Lo and Mei f4(|. Dysthe fT^| pointed out 
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Fig. 7. Evolution of one-soliton solution (15.11) to NLSE (11.11) at kh = 3. SSF parameters: 
At = 1, Af = 2, £ E [-1000, 1000); 1/ « 0.0566 



Fig. 8. Evolution of one-soliton waveform 115. ID taken as the initial condition in HONLSE 
J4.16D at kh = 3. SSF parameters: At = 0.5, A£ = 2, £ E [—4000, 4000). Accuracy: 
| T —50000 : Er(/o) = 0.050% 



£ 


Fig. 9. Variations in the amplitude of the quasi-soliton solu¬ 
tion with distance at kh = 3 



Fig. 10. Mean wave speed as a function of distance at kh = 3: 
solid curve — quasi-soliton, dashed line — NLSE soliton (V ~ 
« 0.0566) 
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Fig. 11. Evolution of one-soliton waveform (15.11) taken as the initial condition in 
HONLSE HI. 161 1 at kh = 2. SSF parameters: At = 0.5, A£ = 2, £ € [—4000, 4000). 
Accuracy: |t=60000 : Er(Jo) = 0.028% 



Fig. 12. Variations in the amplitude of the quasi-soliton so¬ 
lution with distance at kh = 2 


that this phenomenon originates due to the wave- 
induced mean flow, whose component in the direc¬ 
tion of propagation of the wave causes a local Doppler 
shift. Here, we proved for the first time that this well- 
known phenomenon can be observed on finite depth 
as well. This result is the main practical achievement 
of our study. Figures [TT1 |T21 and Il3l demonstrate that 
the same quasi-soliton solution and frequency down¬ 
shift are observed at a smaller depth, kh = 2. 

Finally, the free surface profile reconstructed with 
formula (14.191) is shown in Fig.[l3]for kh = 3. The di¬ 
mensionless maximum free surface elevation is about 
0.046. The case kh = 3 corresponds to wavelengths 
twice as large as depth, A ~ 2 h. The typical depth 
of the shelf near the north-west shore of the Black 
Sea varies from 10 to 100 m. Hence, the wavelength 
corresponding to kh = 3 falls within the range from 



Fig. 13. Mean wave speed as a function of distance at kh = 2: 
solid curve - quasi-soliton, dashed line - NLSE soliton (V ~ 
» 0.0469) 


T= 10000 



Fig. 14- Free surface profile with envelope in the form of 
quasi-soliton at kh = 3 
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Fig. 15. Effect of correction A on the coefficients ai o, 0 and 
« 0 , 0 , 1 



Fig. 16. The left- and right-hand-sides of Eq. (lB2l versus kh 

20 to 200 m, which is quite typical of water waves 
observed on the Black Sea. For h = 30 m, we have 
A ~ 60 m and k = 0.1 m -1 . The corresponding 
maximum free surface elevation of the wave shown in 
Fig. HU is about 0.5 m, and the significant wavetrain 
width is about 2 km. Thus, quasi-soliton solutions ob¬ 
tained in this study can describe swells propagating 
on the relatively calm background on seas with inter¬ 
mediate depths. The typical trough-to-crest height of 
such swells is about 1 m. 

6. Conclusions 

The HONLSE derived earlier by Sedletsky (47:] for 
the first-harmonic envelope of slowly modulated grav¬ 
ity waves on the surface of finite-depth irrotational, 
inviscid, and incompressible fluid with flat bottom 
was rewritten in the dimensionless form suitable 
for numerical simulations. One-soliton solutions to 
NLSE are transformed into quasi-soliton solutions 
with slowly varying amplitude when the HONLSE 
terms are taken into consideration. These quasi- 
solitons represent the secondary modulations of grav¬ 
ity waves. They propagate with nearly constant speed 


and possess the unique property of solitons to ex¬ 
ist over long periods of time without breaking. Their 
speed was found to be higher than the speed of the 
NLSE solitons taken as initial conditions in compu¬ 
tations. This phenomenon was observed earlier both 
in experiment and numerical modeling in the case of 
deep-water limit |S|. It is related to the frequency 
downshift originating due to the wave-induced mean 
flow (l9l . [40| . The quasi-soliton solutions obtained in 
this study describe swells propagating on the rela¬ 
tively calm background on seas with intermediate wa¬ 
ter depth. 

The authors are grateful to Dr. S. S. Rozhkov for 
initial discussions that motivated us to undertake this 
study. D. Dutykh would like to acknowledge the hospi¬ 
tality of Institut fur Analysis, Johannes Kepler Uni- 
versitat Linz, where this work was performed. 

APPENDIX A 

On the Correction Introduced by Slunyaev 

Slunyaev re-derived HONLSE (14.ID and introduced a cor¬ 
rection, 

-a(a 2 - l ) 2 (5cr 4 - 18a 2 -3 )k 2 h 2 + 

+ a 2 (a 2 — l ) 2 (a 2 — 9) kh + a 3 (a 2 — l) (a 2 — 5)^ , (Al) 

to the coefficients Q 41 = ai, 0 , 0 and Q 42 = ao, 0,1 derived ear¬ 
lier by Sedletsky 47]. Actually, this correction was deliberately 
ignored by Sedletsky in view of its smallness. Indeed, Fig. 1151 
proves that A can frankly be ignored at kh > 1. 

APPENDIX B 

Relationship to the Sasa—Satsuma Equation 

Taking into account that (|u| 2 )^ = u^u* +uu|, Eq. J4.16D can 
be rewritten in another form: 

u T = —id 2 U££ + iao, 0 , o\u\ 2 u + 

+ ~ ai, 0 ,oM 2 ^f — n 0 , 0 , in (M 2 )f), (Bl) 

ni, 0,0 = o>i, 0,0 — no, 0 , 1 - 

When 

3(—a3)ao,o,0 = -ai,o,o> (B2) 

Eq. dBlD is reduced to the Sasa-Satsuma equation [2^| , which 
possesses an infinite number of integrals of motion and admits 
some additional exact multi-soliton solutions in contrast to 
HONLSE with arbitrary coefficients 0. However, it is clearly 
shown in Fig. |16l that the above relationship among the param¬ 
eters is not satisfied for any kh. Therefore, the Sasa-Satsuma 
equation cannot be obtained from Eq. H4.16D . 
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APPENDIX C 

Conservation of the Wave Action 

Multiply Eq. dBID by u* and the conjugate equation by u, 
Ut = — ia 2 U££ + iao, o, 0 \u\ 2 u + 

+ — ai, o : o|«| 2 «f — ao, 0,1“ (|«| 2 )f), | X «*, 

u * = ia 2 «| ? — iao, o, ol u | 2 “* + 

+ - “1,0,0 l«| 2 «| - a 0, 0, 1U* (| u| 2 )e), | x U, 

and add these two equations: 

(lL*U T + UU*) = —i-(L2 + 

+ <13 + uu^J — 

— at, 0,0 (\u\ 2 u*U£ + |u| 2 im|) - 2a 0 , o, 1 M 2 (M 2 )? ■ 

After some algebraic transformations, we have 

(M 2 ) t = -i«2 + 

+ “3 - (u 5 u|)^ + (u^)J - 

- - (ai, 0,0 + 2ao, o, l) (M 4 )f. 

In the last term we took into account the following relation 
M 2 (M 2 )e = ««*(««*)? = i(uu*uu*) € = ~(M 4 ) ? . 
Integrating this equation over £ from — oo to oo yields 


oo 

/ (M 2 K 


d£ = 0 O- Iq = 


oo 

/ ltt,: 


df = const, 


(Cl) 


where we used the fact that the function u vanishes at d=oo 
along with its derivatives. 


APPENDIX D 

Split-Step Fourier Technique 

1. Linear equation 

Consider the linear part of HONLSE (14.161 1: 
u r = —ia 2 U££ + a 3 tt£££, u = u(£, r). 

Apply the Fourier transform to the function w(£, r): 


oo 

«(«, r ) = / “(£’ r ) exp(-ire£)d£ = T«[ti(f, r)]. 

— OO 

The inverse Fourier transform is written as 

oo 

ti(£, r) = J u(k, t) exp(iK^)dK = ^ _ 1 [G(«, r)]. 


(Dl) 


(D2) 


(D3) 


The Fourier transforms of the derivatives of function u(£, r) 
are expressed as 

(u^) = i ku, (u^) = —k 2 u, ..., (u n ^) = (i«) n t£. (D4) 

/SSTV 2071-0186 . Ukr. J. Phys. 2011 Vol 59, No. 12 


Hence, linear equation m takes the following form in the 
Fourier space: 

u T = (—ia 2 (itt) 2 + a 3 (itt) 3 ) u, u( 0) = uo. (D5) 

This ordinary differential equation can easily be integrated, 
u = uo exp((ia 2 /^ 2 — ±asK 3 )r) , (D6) 

and the following solution for u(£, r) is obtained: 

oo 

u= J uo exp((ia 2 /“C 2 — ia 3 « 3 )r) exp(i/c£)d/c. (D7) 

— oo 

2. Nonlinear equation 

Nonlinear equation (14.16D can be split into the linear and non¬ 
linear parts: 

u r = —ia 2 U££ + iao, 0 ,Ou|ifc| 2 + 

+ (asu£££ — ai, o, ou^\u\ 2 — ao, o, l u 2 u^j = (C + J\f) u , 
where 

£ = — ia2^ + asdg ££, (D8) 

A/" = iao, o, oM 2 — ai, o, o u £ u * — ao, o, iuug (D9) 

are the linear and nonlinear operators, respectively. The semi¬ 
discretization in time is performed as follows: 

«(£, T + At) - u(£, r) I 

--- = (£ + N) u(£, t) => 

At I At—^0 

=> «(£, t + At) S3 it(£, r) + Ar(£ + Af)u(£, r) ~ 

« e AT(£+A °u(£, r), 

and the second-order Strang formula for noncommuting oper¬ 
ators [Hill is used: 

e Ar(£+A0 = 5( 2 ) (At) = 


= exp (—A/) exp (A tC) exp (-^-A/J 
/At 


(D10) 

(Dll) 


= exp (rfc) exp (ArAf) exp 

In our computations, splitting (|D11|> proved to be more accu¬ 
rate than iIdToI The linear part is integrated exactly using 
relation m 

e Ar£ u(£, t) = 

= T/ 1 [e AT (- 1 ° 2 ( lre ) 2 +' I 3( lK ) 3 )^[u(^, t)]] , (D12) 

and the nonlinear part is corrected at each step as follows: 
e ArJV u(£, t) = 

_ e A T( i ao,°,ol“l 2 - Ili , o ,o“c u *- [lo ,o, iuu i) T ) (D13) 

Following Yoshida (65| . a more accurate fourth-order splitting 
can be introduced as well: 

S (4) (At) = S (2) (piAt) S (2) (p 0 Ar) S (2) (piAr), (D14) 
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V 0 = -7 


2 1 /3 


—1.70, p 1 = 


1 


1.35. 


2 - 2 1 / 3 ’ 2 - 2 1 / 3 
For a more detailed description of the SSF technique, the reader 
can refer to 14311. 
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HEJIIHIHHE PIBHHHHH TTT PE TT THIEPA 

Biimoro nopaflKy jyin obbi^hoi noBuibHO 

MO^yjIBOBAHHX rPABITAU,IHHHX XBHJIb 


HA nOBEPXHI PIflHHH CKIHHEHHOI rjIHBHHH 
TA HOrO KBA3ICO,niTOHHI P03B’H3KH 

P e 3 io m e 

Po3rji5myro HejiimHHe piBHHHHH Ulpe r zj;iHr'epa BHiupro nopH,zj,- 
Ky, BHBe^eHe paHime IO.B. Ce r n,jien;BKHM [YO>K 48(1), 82 
(2003)] ,h,jih obBi^HOi nepnioi" rapMOHiKH noBijiBHO MO^yjibOBa- 
hhx rpaBiTan,iHHHx xbhjib Ha noBepxHi 6e3BHxpoBo'i, HeB’H3KOi 
Ta HeCTHCJIHBOl pi^HHH 3i CKmueHHOIO rjinbHHOIO i njIOCKHM 
flHOM. Up piBHHHHH BpaxoBye ,u,Hcnepciio TpeTboro nopsmicy i 
KybiuHi HejiiHiHHO-,u,HcnepciHHi flO^aHKH. B ^amil poboTi boho 
npHBe^eHO #0 5e3po3MipHoro BHrjm^y, b HKOMy <J>irypye jih- 
me oflHH 5e3po3MipHHH napaMeTp kh, ^e k — xBHJibOBe hhcjio 
Hecyuoi’ xBHjii, a h — He36ypeHa rjinbHHa pi^HHH. IIoKa3aHO, 
iu,o npH BpaxyBaHHi ^o^aHKiB BHiupro nopn^Ky o^HOCOjiiTOH- 
Hi po3b’h3kh KJiacnuHoro HejimmHoro piBHHHHH lUpe^iHr'epa 
nepeTBOpiOIOTbCH B KBa3icOJliTOHHi P03 b’h 3KH 3 nOBijIbHO 3MiH- 
hok) aMnjiiTya,oio. LX,i KBa3icojiiTOHHi po3b’h3kh npeppTaBjiH- 
iOTb BTopHHHi MO,zi,yjiHij,ii rpaBiTau,iHHHx xbhjib. 
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